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Abstract 

We derive a variational formula for the outward normal derivative of the Green func- 
tion for the Schrodinger and Laplace-Beltrami operators, viewed as perturbations of the 
Laplacian. As an application we begin to characterize elliptic growth — the growth of a 
domain pushed outward by its own Green function. 

1 Introduction 

The outward normal derivative of the Green function of a domain is of central importance 
to boundary-value problems associated to the underlying elliptic operator. Given a domain 
D C C, a point w € D, and an elliptic differential operator L, we define the Green function 
g w : D — {w} — > R as the solution (if it exists) to 

{Lg w = S w in D 
g w =0 on dD, 

where S w denotes the Dirac delta supported at w. Occasionally we will denote g w (z) by g(z, w), 
emphasizing the functional dependence upon both inputs. The existence of a Green function 
for a domain is tied to the regularity of the boundary. We will work only with smooth and 
analytic boundaries, which guarantees existence of Green functions and allows use of Green's 
theorem. 

The Green function is used to solve the inhomogeneous boundary-value problem for L with 
zero boundary data. That is, a solution to 




is given by 

u(w) = / g w fdA. 



Of equal importance is the homogeneous boundary-value problem 

{Lu = in D 
1 u = f on dD. 

In some situations, the Green function can be used to extract a solution to this problem as 
well. In the case of the Laplacian, L — A, a solution to ((TJ is given by 

u(w) = / fd n g w ds, 

JdD 
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where d n denotes the outward normal derivative. The integral kernel Pq{w) = d n g w {C) is called 
the Poisson kernel of the region D, and for simplicity of notation we sometimes use P^{w) in 
place of d n g w (£), especially when we wish to view £ as a fixed parameter. 

If instead we consider an operator in divergence form, L = VAV, then a solution to ([1]) is given 

by 

u(w) = / f\d n g w ds. 

JdD 

The normal derivative of the Green function also appears in the dynamics of moving boundaries. 
In [T] the authors define a process wherein a domain grows with a boundary velocity determined 
by the outward derivative of its own Green function. This so-called elliptic growth is used to 
model real processes such as electrodeposition, crystal formation, and quasi-static fluid flowing 
through a porous medium. When the underlying operator is the Laplacian the growth is com- 
monly called Laplacian growth or Hele-Shaw flow — see [4] for a survey. When the underlying 
operator is more general however, much less work has been done. It is known for example that 
lemniscates cannot survive a Laplacian growth process (see [2]), but a carefully-chosen elliptic 
growth might preserve them. To approach such problems we need an understanding of how 
the Green function depends upon the operator. 

Compared to that of a general elliptic operator, the Green function of the Laplacian is well- 
understood. In [3] the author develops formulas for Green functions under various perturba- 
tions of the Laplacian; here we follow the same paradigm, studying the normal derivative of 
the Green function under a perturbation of the operator. 



2 Variation of the Normal Derivative 
2.1 Schrodinger 

Our starting point is the variational formula derived in [3 . Consider a bounded domain D 
with smooth, analytic boundary, a positive function u £ C°°(D) and z,w G D. For e > the 
Green function g^ of A — eu satisfies 

g* w (z) = g w (z) + e f u(£)g(z, £)g(w, dA(0 + o(e) (2) 

J D 

as e — > 0, where g w denotes the Green function for the Laplacian and the error term converges 
uniformly in z for each fixed w. Furthermore, a full series expansion is given by 

oo 

S ; = ^(TM) n S „, 

71=0 

where we've used the operators M : (p i— > wp and 

T : (p *—> / g w <p dA. 
Jd ' 

To derive a similar formula for d n g^ our major tool is the following lemma. 

Lemma 1. Let D C C be a bounded domain with smooth, analytic boundary and f G C 1 (D) f] 
C 2 {D). Suppose further that f = on dD. Then for £ £ dD we have 

^(0 = J D Af(z)P(z,C)dA(z), 

where P(z,C) = d n g z (£) is the Poisson kernel of D. 
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Proof. Let A C dD be measurable with respect to arc length on dD. We define an auxiliary 
function u which solves the Dirichlet problem 

jAu = inD 
1 u = xa on dD 

where \A denotes the characteristic function of A. From the Green identity 

df du s 



L{ u %- f0 £) ds =L iuA, - f * u)dA 



we use the assumptions on u and / to hnd 

df 



dn 



ds 



uAfdA. 



(3) 



(4) 



Given z £ D we return to the Green identity Q with the functions g z and u. It follows that 

dg z 



u(z) 



u— — as = 

dD dn 



dn 



ds. 



Inserting this into (U]) and swapping the order of integration gives 



an 



D 



Since A is arbitrary, we conclude 

df 
dn 



(0 



D 



Af(z)jjL(z,C)dA(z) 



Af(zAzX)dA(z) 
on 



ds(C). 



almost everywhere with respect to arc length. Our assumptions on / imply that each side of 
this equation is a continuous function of £. Hence the equation in fact holds everywhere. □ 



Before we can proceed with to the main result, we need to ascertain the regularity of the 
relevant Green functions. To use lemmaQ]we need to be sure the Green function for Schrodinger 
operators are continuously differentiable in some neighborhood of the closure of the domain. 

Lemma 2. Suppose that D C C is a bounded domain with smooth, analytic boundary. Then 
the Green function of the Laplacian on D is analytic in a neighborhood of D except at the 
singularity. 

Proof. Fix w € D. Since dD is smooth analytic, for each £ e dD there is a neighborhood 
of ( and a conformal map ip : N$ — > B so that <f(N^ D D) = D + , the upper half-disk. Taking 
the union of all such Nq and D gives D, a neighborhood of D. We can extend g w to D as 
follows. 

Given £ € dD find the corresponding neighborhood and conformal map ip. Since g w o ip^ 1 
is harmonic on the upper-half disk, on the interval (—1, 1), and continuous up to the real 
axis, the reflection principle for harmonic functions allows us to extend g w o ip -1 to a harmonic 
function on D via the definition 

9w ° = ~9w ° < P _1 (^) 

whenever z € D is below the real axis. Then g w = (g w oip~ 1 )otp extends to a harmonic function 
on D U iV^. Now suppose that z G C — D is a point in both and N^; we need to show that 
the extensions arising from £ and £ agree at z. Consider the conformal maps (p : — > D and 
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ip : — » D. Then ip o ip 1 is a map preserving the unit disk and the real line; hence for some 
real a E (—1, 1) we can define M a : D — > D as 

t — a 

M a (t) = 



I -at 

so that ip o ip~ x = M a . We can rewrite the extension of g w resulting from as 



{g W O = (S« ° ¥> ° -^a~ )(-^a(v( z ))) 



= ( 5w o^- 1 oM- 1 )(M a (^(2))) 
= {g w ° ^ _1 )(</?(z)), 

as desired. Since f/ w is harmonic away from the singularity, it is analytic there. □ 

Lemma 3. Suppose that D C C is a bounded domain with smooth, analytic boundary. Let u 
be a positive, smooth function in a neighborhood of D. Then the Green function of D for the 
operator A — u is smooth in a neighborhood of D except at the singularity. In particular, the 
Green function satisfies the conditions of lemma [7J 

Proof. Fix a point w G D. We can define 

T : ip ^ / g z <pdA 

JD 

which as a left inverse to the Laplacian is a map from the Sobolev space H k (D) to H k+2 (D) 
for any k > 0. Define M to be the multiplication operator arising from u. Since u is smooth 
in a neighborhood of D, the operator M maps H k (D) into itself. It was shown in [3] that 

oo 

g* w = Y J {TMTg w (5) 

71 = 

holds throughout D. Since g w is smooth in D — {w}, so is g^. Furthermore, since g w extends 
smoothly into a neighborhood of D excluding w, equation ([5]) gives a smooth extension of <?* 
into the same domain. □ 

Theorem 4. Let D C C be a bounded domain with smooth, analytic boundary and fix w £ D. 
Suppose that u € C°° (D) is a positive function. The outward normal derivative of the Green 
function g* of the Schrddinger operator A — eu satisfies 

dngUC) = d n gw(0 + e / u(z)g(z, w)P(z, alA(z) + o(e) 

JD 

as e — > 0, where the convergence of o(e) is uniform in £ /or eac/i /ised iu. 
Proof. Notice that g^, satisfies 

(Aff™ = o"™ +e"ff™ in.D 
\g* = on dD 

Using lemma [T] gives 

d n g* w (0 = f (S w + eu(z)g* w {z)) P(z, Q dA(z) 

JD 

= d n g(C,w) + e[ u(z)g* w (z)P(z,C)dA(z), (6) 
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where we have used the fact that P(z, — d n g z (0- Inserting the perturbation formula of 
equation ([2]) into ([6]) gives 

d n g* w (0 = d n g((,w)+e f u(z)g w (z)P(z,() dA + o(e). □ 

We remark that this formula can be easily derived from formal manipulation; beginning with 
equation ([2]) take the outward normal derivative in the z variable at a point £ on the boundary. 
Assuming o(e) remains small and the derivative commutes with the integral, we produce the 
result. The above proof actually produces an expression for the second variation of <9„g* u , 
which we give here as a corollary. 

Corollary 5. Under the same conditions as the previous theorem, the second variation ofd n g^ 
is given by 

S 2 d n gZ(0= ! I u(z)u(0g(t,w)g(t,z)P(z,0dA(0dA(z). 

J D J D 

As a further corollary, we can deduce a monotonicity result. 

Corollary 6. Under the same conditions as the previous theorem, the first variation 5d n g^(^) 
is nonpositive at each point of dD. 

Proof. For any z,w £ D and £ G dD we have P(z, C), u{z) > and g(z, w) < 0. Thus 

/ u(z)g(z, w)P{z, C) dA(z) < 0. 
Jd 

The result follows from theorem [4] □ 
2.2 Laplace— Beltrami 

Having completed the analysis of the Schrodinger case, we turn to a different perturbation of the 
Laplacian. As before, consider a bounded domain DcC with smooth, analytic boundary and 
u G C°°(D). Given e > define A = 1 + eu and the Laplace-Beltrami operator L = VAV. Let 
g^j denote the Green function for L with singularity at w € D; we wish to derive a perturbation 
formula for (? n <7* (£) when e is close to 0. The result we obtained for the Schrodinger operator 
turned out to be the obvious one that results from formal manipulation. In this vein, we can 
anticipate a Laplace-Beltrami analogue by using the following perturbation formula from [3]: 

gl( z ) = gw(z) + e uVg z -Vg w dA + o(e) (7) 
Jd 

as e — > 0, where w, z € D and o(e) converges uniformly in z for each fixed w. Let's take a 
normal derivative in z at the point £ G dD, assuming that o(e) remains small and derivative 
commutes with integral. We are lead to the formula 

dngtiC) = d n g w (C) + e / «VP ( • Vg w dA + o(e). 
Jd 

Unfortunately this integral might diverge, but formal manipulation and integration by parts 
leads us to the following theorem. 

Theorem 7. Let D C C be a bounded domain with smooth, analytic boundary and fix w G D. 
Suppose that u G C°°{D) is a positive function and define X(z) = 1 + eu(z) for e > 0. The 



5 



outward normal derivative of the Green function g* of the Laplace- Beltrami operator L — VAV 
satisfies 



d n 9w(0 = 9„g w (0 + 



Aug w PQ dA - d n g w (() [u(() + u(w)] 



D 



o(e) 



as e — > 0, where the convergence of o(e) is uniform in £ for each fixed w. 



Proof. For a rigorous proof, we use a well-known change of variables to reduce the problem 
to that of the Schrodinger operator. Since u € C°°(D) C L°°(D), we consider only those 
e satisfying < e < eo < HuH^ 1 . This ensures that A > and that all of the following 
linear approximations have uniformly small error terms. If we define the auxiliary functions 
G*„ = g%\X and 

aVa 



v = 



Va 



then we claim •v/ X(w)G^ is the Green function for A — V . Verification is a straight forward 
computation which we omit; a full derivation is given in (3J. Furthermore, the dependence of 
V upon e can be seen to be V — eAu/2 + o(e). Using theorem 2] we obtain 



= d n g w (() + ~ J Aug w P c dA + o(e). 



^/\(z)\(w)g* w (z) 
Using g^(() — the left-hand side of this equation simplifies to 



(8) 



y / X(z)\(w)g* w (z) 



9*w(C)d„ 



V / A(z)A(u-) 



^\(()\(w)d n g* w (C) 



d n 9w(0 Au 9 wP C dA + °( e ) 



- ^\(C)\(w)d n gU0 

We can rewrite (j8]) as 

d n g* w (0 = WC)AH)" 1/2 
We only want first-order dependence upon e, so we can expand 

(»))^ = l-e(^) +0 (e) 

so that (0) becomes 

'u(C) + u(w) 



(9) 



1 -e 



o(e) 



dn9w(() + 2 / Au 9™ P C dA + °( e ) 



9n9w(0 + 



D 



Aug w P c dA - d n g w (0 l u (0 + u ( w )} 



as desired. 



□ 



As we did for the Schrodinger operator, we can deduce another monotonicity result. 

Corollary 8. Under the same conditions as the previous theorem, the variation Sd n g^,(^) 
satisfies 

5dn9* w ds < 0. 

dD 
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Proof. From the previous theorem, 



5d n glj ds — Aug w dA — u(w) — / ud n g w ds, 

dD JD JdD 

where we have swapped integrals and used the fact that 

/ d n g w ds= / P ( (z)ds(() = l 

JdD JdD 

for each w, z € D. The integral swap is valid since all integrands are of constant sign. Green's 
identity allows us to write 



/ 5d n g* w ds = -2 / ud n g w ds < 0, 

JdD JdD 



IdD JdD 

since u,d n g w > everywhere on 3D. □ 

3 Applications 

3.1 The Dirichlet Problem 

For an elementary application of the Green variation, consider solving the Dirichlet problem 
on a domain D: 

((A-eu)ip e =0 inD 
| (f e = f ondD 

where / is a bounded Borel function, u > is smooth in a neighborhood of D and e > is 
small. If we ignored the eu term in the problem — that is, assume e is zero — we can estimate 
the error with the following theorem and corollary. 

Theorem 9. Let D C C be a bounded domain with smooth, analytic boundary and suppose 
u G C°°(D) is a positive function. For e > and a bounded Borel function f the solution ip £ 
to the Dirichlet problem (|10[) satisfies 



(p e (z) = (p (z) + e I wp g z dA + o(e) 

J D 

as e — > 0, where the error term converges uniformly in z. 

Proof. Let g* u denote the Green function of the operator A — eu in the region D. Note that 

J (A - eu)((p e ~ ifo) = tutpo inD 
1 ip e — ipo = on dD, 

whence an expression for tp € is given by 

<Pe(z) = Vo(z) + e / u(p g*dA. 

J D 

From the perturbation formula ([5]) we have 

ife(z) = (fio(z) +e uipogz dA + o(e), 
J D 

as desired. □ 
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Corollary 10. With the same assumptions as the previous theorem, the linearization of (p e 
has the pointwise bound 

\6<p e (z)\ < NI2IM2II/H00. 

Proof. From the maximum modulus principle for harmonic functions, sup D \<po\ < sup aD |/|. 
The result follows from this and the Cauchy-Schwarz inequality. □ 

3.2 Elliptic Growth and an Inverse Problem 

Consider a bounded domain DcC with smooth, analytic boundary and fix w £ D. Given a 
positive function u € C°°(D), an elliptic growth process governed by the Schrddinger operator 
A — u is a nested sequence of domains {D(t) : t > 0} such that D(0) = D and 

V n (0 = d n g* w ((), 

where V n denotes the outward normal velocity of the boundary dD(t) and g* w denotes the Green 
function of the operator A — u for the domain D(t). That is, the Green function "pushes" the 
boundary of the domain outward, changing the domain — and hence the Green function — in a 
nonlinear, nonlocal manner. In beginning to understand the nuances of this process, we ask the 
following question. Which families of nested, growing domains can be produced by an elliptic 
growth process? 

Our first approach to addressing this question is to consider the local problem. Given the initial 
domain D we seek which velocity fields V can emerge from different choices of the function u. 
The intuition is as follows. Let Vb be the velocity field that arises from elliptic growth with 
u = (this is the well-known Laplacian growth process). If V is another velocity field close to 
Vo, can we choose u ~ so that V arises from elliptic growth via A — u! In other words, if we 
view the variation as a linear map u i-> Sd n g£, on an appropriate function space, what is the 
image of this map? We are led to the following theorem. 

Theorem 11. Let D C C be a bounded domain with smooth, analytic boundary and fix w e D. 
Define 

Au(Q= f u(z)g(z,w)P(z,QdA(z). 
Jd 

Then A is a linear map L 2 (D) — > L 2 {dD) with dense range. 

Proof. To prove A is defined on all of L 2 (D) it suffices to show that g w P(; € L 2 (D) for each 
w e D and £ e 3D. Away from w and ( the function g w P( is smooth; furthermore, the 
singularity at w is logarithmic, which is square-integrable. It remains to understand the 
singularity at £. Having assume the boundary of D is smooth and analytic, there is a conformal 
map from a neighborhood of £ to the unit disk which transforms the Poisson kernel and Green 
function canonically. Hence it suffices to understand the singularity at £ in the case when 
D = D, the unit disk. Without loss of generality, we can assume w = temporarily; a Mdbius 
transform can preserve the disk and send w to 0, which doesn't disturb the behavior of the 
singularity at £. Note that 

90(Z) P M ^(<±£). 

Taking the limit z — > £ gives 



lim g (z)P(z,Q = - 




so goP^ remains bounded near £. At this point we return to the general domain D; we can 
write g w (z)P(z, £) = h^(z), a bounded function in a neighborhood of D except at w, where it 
has a logarithmic singularity. We conclude g w P^ € L 2 (D). 
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Next we show that A maps into L 2 (dD). Note that 



/ \Au\ 2 ds< f \\u\\ 2 \\h c \\ 2 ds(C) < \\u\\ 2 sup \\h c \\ 2 ■ length(0£>), 

JdD JdD C&9D 

so it suffices to show that sup aD \\h^\\ 2 < oo. This follows easily from the fact that £ i-> ||/i^||2 
is a continuous function on the compact set <9£>. 

Next we show that the range of A is dense in L 2 (dD). Suppose v € L 2 (dD) is chosen so that 

/ v (C)(A«)(C)d*(C)=o 

for all u G L 2 (D). Using the definition of A and swapping integrals gives 

J u{z) Qf w(C)M*) <MC)) <*A(*) = 0. 

Since it is arbitrary, the inner integral must vanish for almost all z € D. Note that = 
g w (z)P(z, Q and ^(z) < for all z G D, so 

/ «(C)P(«,C)d*(C)=0 

for almost all z E D. This integral is a harmonic function in z which is zero almost everywhere 
in D, hence everywhere by continuity. The domain is regular enough so that, when z approaches 
dD, we recover the function v. Hence v — in L 2 (dD). This shows the orthogonal complement 
of the range of A is trivial, so the range of A is dense in L 2 (dD). □ 

For physical reasons, we might prefer function spaces different from L 2 . One natural candidate 
for velocity fields is L 1 (9_D), since the L 1 distance between velocity functions measures the 
difference in areas due to the domain growing in different ways. Furthermore we could consider 
only smooth perturbation functions u, in the spirit of the original perturbation result. To these 
ends we have the following corollary. 

Corollary 12. With the same hypotheses as the previous theorem, A maps C°°(D) — > L 1 (9_D) 
and has dense range. 

Proof. Since dD has finite measure, L 2 {dD) is dense in i 1 (9_D). Furthermore C°°(D) is dense 
in L 2 (D), and the result follows. □ 

The inverse problem is not fully solved with these results. First we should only consider 
nonnegative functions u and characterize the range of the map A. Secondly, the restriction of 
the domain of A to smooth functions begs the following question. Can every smooth function 
in i 1 (9Z?) be obtained from A when restricted to C°°(D)1 Finally, which of these results can 
be carried over to elliptic growth governed by a Laplace-Beltrami operator? In that setting the 
relevant integrals are more delicate and require further analysis. All of these problems provide 
directions for future work. 
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